Correlation induced half-metallicity in a ferromagnetic single-layered compound: 

Sr 2 Co0 4 



O 

(N 

G 

>—> 

On 



Sudhir K. PandejQ 

UGC-DAE Consortium for Scientific Research, University Campus, 

(Dated: January 9, 2010) 



Khandwa Road, Indore - 452001, India 



The electronic and magnetic properties of S^CoCU compound have been studied using ab initio 
electronic structure calculations. As opposed to GGA calculation, which gives ferromagnetic metallic 
solution, GGA+C/ calculations provide two kind of ferromagnetic solutions: (i) half-metallic and (ii) 
metallic. The half-metallic solution is a ground state of the system and the metallic one is a 
metastable state. The strong hybridization between Co 3d and O 2p orbitals decides the electronic 
and magnetic properties of the compound. The total magnetic moment per formula unit is found to 
be ~ 3 fiB (S = 3/2). Our calculations give the magnetocrystalline anisotropy energy of ~ 2.7 meV, 
which provides a good description of experimentally observed large magnetocrystalline anisotropy. 
The Heisenberg exchange parameters up to fourth nearest neighbours are also calculated. The 
mean-field theory gives the Tc = 887 K. The possible physical implications of the ferromagnetic 
half-metallic ground state are also discussed. 

PACS numbers: 71.27.+a, 71.20.-b, 75.20.Hr 
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I. INTRODUCTION 

The single-layered compounds with a general formula 
of A2BO4 (A and B stand for rare-earth/alkali metals and 
transition metals, respectively) have attracted a great 
deal of attention after the discovery of high temperature 
superconductor in I^-ajSr^CuC^.— These compounds 
show many exotic physical phenomena like spin/charge 
stripes formation in nickelates and manganitesi; and 
spin-triplet superconductivity in ruthenates.— 

Among the single-layered compounds recently synthe- 
sized Sr 2 Co04 has been reported to show ferromagnetic 
(FM) and metallic behaviours, which have never been 
found in any other such materials4^ The magnetiza- 
tion data show Tc around 250 K. The large magnetic 
anisotropy where magnetic easy axis is the c-axis is also 
observed. Moreover, the resistivity data reveal the quasi- 
two-dimensional electronic nature for the system. In- 
terestingly, there are contradicting reports on the mag- 
netic moments. The Matsuno et al£ and Wang et alA 
have reported the saturation magnetization of about 1.8 
and 1.0 /ib/Co, respectively. The effective magnetic mo- 
ment of 3.72 \xb has also been found by fitting Curie- 
Wiess behaviour in the paramagnetic phase, which sug- 
gests the magnetic moment corresponds to S ~ 3/2 spin 
configuration^ 

Lee et al^ tried to understand the above contradic- 
tory reports on the saturation magnetization by using 
LDA+[/ calculations. Their studies suggest that such 
values of magnetization can be intrinsic to the system as 
they got similar two values for the total magnetic mo- 
ment for different ranges of U . However, if we read the 
works of Matsuno et al. and Wang et al. carefully we 
can easily make out that the different values of satura- 
tion magnetization they got may not be intrinsic to the 
system. The M(H) data of both the groups look similar. 
For example, the values of total magnetization at 5 T 
found by Wang et al. and Matsuno et al. are about 1.4 



and 1.6 /ib/Co, respectively, which are not very different. 
The different values of saturation magnetization quoted 
in the last paragraph appear to arise due to two different 
approaches the authors used in extracting them. 



The partial replacement of Sr by rare-earth elements 
leads to the formation of interesting magnetic phases 
keeping the crystal structure intact^— For example, 
Sri. 5 Lao.5Co04 compound retains the FM ground state 
of Sr2Co04, whereas SrLaCo04 manifests spin glass 
stated Contrary to FM state for Sri^Lao.sCoCU, the 
ground state of Sri.sPro.sCoCU compound is found to be 
spin glass. 7 Similarly, the formation of spin glass state is 
also reported in Sri.25Ndo.75Co04 compound^ The ex- 
istence of spin glass state in the doped compounds ap- 
pears to be a generic phenomenon. It is well known that 
the spin glass state arises due to the presence of com- 
peting ferromagnetic and antiferromagnetic interactions. 
Therefore, it would be interesting to calculate the na- 
ture of magnetic interactions between the neighbouring 
Co atoms in the Sr2Co04 compound. 



Here, we report the detailed electronic and magnetic 
states of Sr2CoC>4 using ab initio electronic structure cal- 
culations. The ground state solution is found to be a fer- 
romagnetic half- metallic state having .5 = 3/2 spin con- 
figuration, which is different from that reported by Lee 
et al. The strong hybridization between Co 3d and O 2p 
orbitals decides the electronic and magnetic properties of 
the compound. The calculation gives magnetocrystalline 
anisotropy energy of ~ 2.7 meV. The nature of mag- 
netic interactions between the Co atoms is calculated up 
to fourth nearest neighbours, which shows mixed ferro- 
magnetic and antiferromagnetic coupling. The estimated 
ferromagnetic transition temperature within mean-field 
theory is ~ 887 K. 
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II. COMPUTATIONAL DETAILS 

The nonmagnetic and ferromagnetic electronic struc- 
ture calculations of Sr2Co04 compound were carried 
out using LmtArt 6.61.— For calculating charge density, 
full-potential linearized Muffin-Tin orbital (FP-LMTO) 
method working in plane wave representation was em- 
ployed. In the calculations, we have used the Muffin- Tin 
radii of 2.915, 1.965, 1.608, and 1.778 a.u. for Sr, Co, 
01, and 02, respectively. The charge density and effec- 
tive potential were expanded in spherical harmonics up 
to I = 6 inside the sphere and in a Fourier series in the 
interstitial region. The initial basis set included 5s, 4p, 
and Ad valence, and 4s semicore orbitals of Sr; 4s, 4p, and 
3d valence, and 3p semicore orbitals of Co, and 2s and 2p 
orbitals of O. The exchange correlation functional of the 
density functional theory was taken after Vosko et alJ£ 
and GGA calculations were performed following Perdew 
et al^- 

The effect of on-site Coulomb interaction (£/) under 
GGA+U formulation of the density functional theory is 
also considered in the calculations. The detailed descrip- 
tion of the GGA+U method implemented in the code 
can be found in Rcf. [12]. The double counting scheme 
used in the code is normally called as fully localized limit 
in the literature. To study the role of orbital degrees of 
freedom, we have also included spin-orbit coupling (SOC) 
in the calculations. The Self-consistency was achieved 
by demanding the convergence of the total energy to be 
smaller than 10~ 5 Ry/cell. (10, 10, 10) divisions of the 
Brillouin zone along three directions for the tetrahedron 
integration were used to calculate the density of states 
(DOS). 

III. RESULTS AND DISCUSSIONS 

The atomic arrangement in the unit cell is shown in 
Fig. 1, which displays bodycenter tetragonal lattice. The 
Co atoms occupy the corners and bodycenter positions 
and each Co atom is surrounded by six O atoms forming 
a distorted octahedron. In this structure there are two 
kinds of O represented by 01 and 02. 01 lies in the ab 
plane and 02 along the c-axis. The bond distance be- 
tween the corner and bodycenter Co atoms is almost 1.8 
times larger than that between nearest Co atoms sitting 
at the corners. This suggests the quasi-two-dimensional 
nature of the system where its electronic properties are 
expected to be decided by transport of electrons in the 
ab plane. 

In order to know the crystal-field effect and nature 
of Co-0 bonding in Sr2Co04 we have plotted the par- 
tial density states (PDOS) of Co 3d and O 2p obtained 
from nonmagnetic solution in Fig. 2. In the octahe- 
dral symmetry Co 3d states split into ti g and e g states 
and the separation between these states is found to be 
~ 1.5 eV. The degeneracies of these states are further 
lifted in tetragonal symmetry as evident from Figs. 2(a) 




FIG. 1: (Color online) Atomic arrangements in the unit cell. 
Sr, Co and O atoms are represented by spheres with decreas- 
ing radii. 



and (b). The triply degenerate t^ g states split into dou- 
bly degenerate (d xz , d yz ) states and nondegenerate d xy 
state. Similarly, doubly degenerate e g states split into 
nondegenerate d x 2_ y 2 and d z 2_ r 2 states. The O 2p or- 
bitals also split in doubly degenerate (p x , p y ) orbitals and 
nondegenerate p z orbital as evident from Figs. 2(c) and 
(d). 

The t2g and e g sectors are spread over the energy range 
of about 5 and 10 eV, respectively. Roughly two times 
larger extent of the e g states is due to larger overlap be- 
tween e g and p orbitals in comparison to that between ti g 
and p orbitals. In the present situation one can get max- 
imum overlaps of d x i_ y i and 01 (p x , p y ) orbitals; and 
d z 2_ r 2 and 02 p z orbitals which can lead to larger band- 
width of e g symmetric states. The occupied d z 2_ r 2 states 
are mostly found in the narrow region of 1 eV whereas 
occupied d x 2_ y 2 states are spread over an energy window 
of about 4 eV. Total number of d electrons is found to be 
~ 6.5, which is about 1.5 more than the expected nom- 
inal value for Co 4+ ion. The above observations can be 
considered as a signature of the covalent nature of the 
Co-0 bonds. There are large Co 3d PDOS at the E F (~ 
8.4 states/eV/atom) in the nonmagnetic solution which 
may be considered as a signature of ferromagnetic ground 
state based on the Stoner theory. 
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FIG. 2: (Color online) The partial density of different 
symmetric states obtained from nonmagnetic solution, (a) 
d X z/dy Z (thin lines) and d xy (thick lines) states, (b) d x 2_ y 2 
(thin lines) and d z 2_ r 2 (thick lines) states. The p x /p y (thin 
lines) and p z (thick lines) symmetric partial density of states 
of planar oxygen (01) and apical oxygen (02) are shown in 
(c) and (d), respectively. 



FIG. 3: (Color online) The density of states obtained from 
ferromagnetic GGA calculation, (a) total density of states 
(TDOS), (b) d X z/dyz (thin lines) and d xy (thick lines) states, 
and (c) d x 2_ y 2 (thin lines) and d z 2_ r 2 (thick lines) states. 
In (d) partial density of p symmetric states for planar oxygen 
(Ol) and apical oxygen (02) are represented by thin and thick 
lines, respectively. 



In order to know the electronic and magnetic struc- 
tures of the compound we have carried out FM calcu- 
lation. The energy of the FM solution is found to be 
~ 336 meV less than that of nonmagnetic solution in- 
dicating the ferromagnetic ground state for the system. 
This result is in accordance with the experimental obser- 
vation and earlier theoretical calculations^ The total 
and partial density of states obtained for the FM config- 
uration are plotted in Fig. 3. There is an asymmetric 
distribution of states in both the spin channels. The dif- 
ference between number of electrons in up and down spin 
channels is a measure of net magnetization and found to 
be ~ 2. Under the rigid band picture one can estimate 
the effective value of the exchange parameter (J) by con- 
sidering the energy difference between the up and down 
band edges, which is estimated to be ~ 0.4 eV. The t 2g 
up-spin bands are fully occupied whereas t 2g down-spin 
and e g bands are partially occupied. 

Similar to the nonmagnetic case, e g bands are highly 
extended in comparison to t 2g bands. The strength of ex- 
change parameters for t 2g and e g electrons are estimated 



to be ~ 0.6 and 0.4 eV, respectively. The larger value 
of J for t 2g electrons can be attributed to more localized 
t 2g states leading to higher Coulomb interaction. Inter- 
estingly, O 2p PDOS also show spin polarization and the 
estimated value of J for p electrons is also found to be ~ 
0.4 eV. The values of magnetic moment for Co, Ol and 
02 atoms are found to be about 1.52, 0.17 and 0.07 /ib, 
respectively. The total magnetization per formula unit is 
~ 2 hb, which is almost same to the earlier reports^ 
The appearance of magnetic moments at the O sites may 
be attributed to the strong hybridization between d x 2_ y 2 
and Ol (p x , p y ) orbitals; and d z 2_ r 2 and 02 p z orbitals. 
The values of magnetic moment for Ol and 02 atoms 
also appear to support this conjecture as the magnetic 
moment of Ol is roughly 2 times larger than that of 02. 

In order to know the role of orbital degrees of free- 
dom on the magnetic properties of the system we have 
included SOC in the calculations by considering magneti- 
zation directions along (001) and (100) axes. Our calcu- 
lations indicate that the c-axis is an easy axis as the total 
energy of the system for (001) direction is found to be ~ 
2.7 meV/Co less than that for (100) direction. This result 
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is in line with the experimentally observed magnetic easy 
axis. 3 This value of magnetocrystalline anisotropy energy 
(MAE) is ~ 2.1 meV/Co more than that obtained by Lee 
et al£- Based on the magnetization data of Ref . [3] , Lee et 
al. have given a rough estimate of ~ 1.5 meV/Co for the 
MAE. This value is about three times larger than that of 
the Lee et al. and about half of our value. The inclusion 
of SOC keeps the spin magnetic moments almost intact 
and induces 0.1 us orbital magnetic moment at the Co 
site. 

The above GGA calculations provide the good repre- 
sentation of experimentally observed electronic and mag- 
netic properties of the system including magnetocrys- 
talline anisotropy. However, the on-site Coulomb inter- 
action parameter U has been found to be important in 
understanding the detailed electronic and magnetic prop- 
erties of the 3d transition metals oxies.— Therefore, it 
would be interesting to see the effect of U on electronic 
and magnetic states of the compound under GGA+J7. In 
this approximation U and J are used as parameters. 

The values of U and J required to reproduce the elec- 
tronic structure of a compound are sensitive to the ap- 
proximations used in the calculations. 13,14 Moreover, it is 
well known that the GGA+/7 calculations often converge 
to the local minima depending on the starting electronic 
configurations usedJ£~— In the light of these facts we 
have varied U from 3 to 5 eV and fixed the value of J 
to 0.4 eV as estimated for the e g electrons. To know the 
exact ground state we have also used different starting 
configurations for the Co 3d electrons. This range of U 
is found to provide the good representation of the exper- 
imentally observed electronic and magnetic properties of 
different Co based systemsi 13 ' 18 " — 

Based on symmetry consideration we used 10 differ- 
ent starting electronic configurations for the Co atom 
in the calculations. For these configurations only two 
kind of FM solutions SI and S2 are found to exist. The 
SI and S2 give half-metallic and metallic states, respec- 
tively. The difference in the electronic structure of Co 
in both the solutions can be found from Fig. 4, where 
we have plotted the Co 3d PDOS calculated for U = 4 
eV. The half-metallicity with a gap of about 1 eV in the 
down-spin channel is evident for the SI. The S2 shows 
metallic state with a large ti g symmetric states around 
the Ep in the down-spin channel. In the up-spin channel 
both the solutions provide similar PDOS of e g character 
in the vicinity of Ep. The comparison of Co 3d PDOS 
obtained from the S2 and FM GGA solutions indicates 
that S2 is originating from the FM GGA solution. The SI 
solution is quite different and it has origin to the orbital 
polarization. 

The PDOS closer to Ep along with occupancies of the 
Co 3d orbitals listed in Table 1 can be used to charac- 
terize both the solutions. Numbers written in normal 
and italics are corresponding to the SI and S2, respec- 
tively. The tig^ states are fully occupied in both the 
solutions. The main difference between the two solutions 
sector. In the SI and S2, d xy ^ state 
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FIG. 4: (Color online) The partial density of Co 3d symmet- 
ric states obtained from ferromagnetic GGA+U (U — 4 eV 
and J = 0.4 eV) calculations corresponding to two conversed 
solutions represented by SI and S2. tig and e g orbitals are 
denoted by solid and dashed lines, respectively. 



can be considered as fully occupied and unoccupied, re- 
spectively. The contribution of d xz ^/d yz ^ state below 
the Ep is small (~ 0.17 electron) and can be considered 
as fully unoccupied in the SI, which gives rise to the 
half-metallic state. However, d xz ^/d yz ^ state is roughly 
half-filled (~ 0.6) in the S2 and mainly responsible for 
the metallic state. Thus, we can represent the SI and S2 
by d X y^ and d X7 } x l 2 / 'dy^ 1 / 2 electronic configurations, 
respectively. 

Table 1: Electronic occupancies of different Co 3d or- 
bitals along with the total number of d electrons for both 
the spin channels obtained from ferromagnetic GGA+U 
(U = 4 eV and J = 0.4 eV) calculations corresponding to 
two self consistent solutions represented by SI (normal) 
and S2 (italics). 



Up 
Down 



d 'XZ I d y Z 



0.95,0.94 
0.17,0.61 



0.94,0.95 
0.92,0. 14 



Up - Down 0.78,0.55 0.02,0.81 0.62,0.65 0.52,0.57 2.72,2.69 



d. r 2_, 



0.84,0.90 
0.22,0.25 



0.88,0.82 
0.36,0.25 



Tot al d 
~4~56,4.55 
1.84,1.86 



arises from the t2 g ^ 



Both the solutions give roughly same magnetic mo- 
ment (~ 2.7 hb) for the Co. The SI induces magnetic 
moment of ~ 0.08 hb at the 02 site whereas the S2 
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provides ~ 0.18 at the 01 site. The total magnetic 
moment per formula unit for both the solutions is found 
to be ~ 3 /zs/fu representative of S = 3/2 spin state. In- 
terestingly, the energy of the SI is ~ 234 meV less than 
that of the S2 indicating that the SI corresponds to the 
true ground state of the system. Present work clearly 
shows the importance of U in stabilizing a half-metallic 
state which turns out to be the lowest energy state of the 
system. Thus, the ground state of Sr 2 Co04 is found to 
be half- metallic with S — 3/2 spin configuration. 

Our half-metallic ground state is quite different from 
that obtained in earlier work£ The half-metallic solution 
of Lee et al. gives S — 1/2 spin configuration. This so- 
lution corresponds to almost occupied d xz ^/d yz ^ states 
and half-filled d xy ^ state whereas our SI corresponds to 
the opposite situation where d xz ^/d yz ^ states are almost 
vacant and d xy ^ state is filled. Here it is worth mention- 
ing that we have also included the generalized gradient 
approximation in the calculations whereas Lee et al. have 
only considered the local density approximation (LDA). 
This may give an impression that the different formula- 
tions of the exchange correlation functional can be the 
reason for the difference between our results and that of 
the Lee et al. To address this issue we have also carried 
out detailed calculations using spin polarized LDA+C7 
method. The results obtained from these calculations 
are found to be almost the same as those obtained by 
GGA+U. This indicates that the existence of two solu- 
tions has nothing to do with the approximations involved 
in treating different exchange correlation functional and 
solely related to the orbital polarization. 

The total and partial densities of states correspond 
to the half- metallic solution are shown in Fig. 5. The 
large spin polarization is evident from the total den- 
sity of states plotted in Fig. 5(a). The lowest occupied 
band shows the signature of van Hove singularity aris- 
ing from the e g symmetric states as evident from Fig. 5 
(c). The d xy orbital is fully occupied and a gap of ~ 1 
eV is created between it and the unoccupied d xz /d yz or- 
bitals in the down-spin channel. In the up-spin channel 
only e g sector provides the conducting electrons with a 
dominating contribution from the d z 2_ r 2 electrons. The 
d x 2_ y 2 symmetric PDOS in the vicinity of the Ep is al- 
most flat suggesting that the temperature dependent be- 
haviour of the compound should be decided by d z 2_ r 2 
electrons. Fig. 5(b) shows a huge exchange splitting (~ 
7.5 eV) for the degenerate d xz /d yz orbitals. This result 
also suggests that the Coulomb correlation in the pres- 
ence of strong hybridization between Co 3d and O 2p 
orbitals along with a tetragonal distortion makes full oc- 
cupation of nondegenerate d xy orbital energetically more 
favorable. This may be considered as a signature of a 
ferrodistortive Jahn- Teller ordering. 

In order to study the detailed nature of magnetic inter- 
actions we have calculated the Heisenberg exchange pa- 
rameters Joi between the Co atoms using FLEUR code^ 
where and i stand for the central atom and coordina- 
tion shell index, respectively. The Co atom sitting at the 
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FIG. 5: (Color online) The density of states corresponds 
to the ferromagnetic half-metallic solution obtained from 
GGA+J7 {U = 4 eV and J = 0.4 eV) calculation, (a) to- 
tal density of states (TDOS), (b) d xz /d yz (thin lines) and d xy 
(thick lines) states, and (c) d x 2_ y 2 (thin lines) and d z 2_ r 2 
(thick lines) states. In (d) partial density of p symmetric 
states for planar oxygen (Ol) and apical oxygen (02) are 
represented by thin and thick lines, respectively. 



origin is considered as the central atom, Fig. 1. The 
values of Joi up to 4 th nearest neighbours (coordination 
shells) are shown in Table 2. It is evident from the table 
that the first and third nearest neighbour interactions are 
ferromagnetic whereas second and fourth nearest neigh- 
bour interactions are antiferromagnetic in nature. The 
value of Joi is found to be f 9.06 meV which is about 15 
times larger than the values of neighbouring interactions. 
This suggests that the effective magnetic interaction in 
Sr2Co04 is ferromagnetic in nature, which is in accor- 
dance with the experimentally observed ferromagnetic 
ground stated Under the mean-field approximation we 
calculated the Curie temperature (Tq) of ~ 887 K, which 
is roughly 3.5 times larger than the experimentally mea- 
sured value. Such discrepancy between the experimental 
and calculated Tq is not surprising as the overestimation 
of Tc under mean-field approximation is well known in 
the literature i22r— 

Table 2: The Heisenberg exchange parameters (Joi) 
up to fourth nearest neighbours (coordination shells) ob- 
tained by ferromagnetic GGA calculations using FLEUR 
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code. and i stand for central atom and coordination 
shell index, respectively. 



Nearest neighbour 


Coordination number 


J i (meV) 




4 


19.06 


2 nd 


4 


-1.11 


3 rd 


8 


1.25 




4 


-0.33 



Finally, we discuss the physical implications of our re- 
sults. The half-metallic ferromagnet has potential appli- 
cation in the sprintonics industry. 2 - This compound can 
be considered as a good candidate as it shows T c closer 
to the room temperature. The d PDOS at the Ep are 
quite low (0.19 states/eV/atom) as evident from the SI 
of Fig. 4. This suggests that even a moderate hole dop- 
ing at the Co sites can lead to insulating ground state. 
The present work predicts the saturation magnetization 
of - 3 /is/fu. The 5K M(H) data of Wang et al. do not 
show any sign of saturation at maximum measured field 
of 5 T where the value of magnetization is found to be 
about 1.4 hb/CoA However, their M(T) data suggest to 
support our calculated S = 3/2 spin state as they have 
reported the effective magnetic moment of ~ 3.72 (xb- 
Here it is important to note that M(H) data of Matsuno 
et al. show almost saturation of magnetization at 7 T 
where the value of magnetization is found to be about 
1.8 /is/Co4 However, this work has been carried out on 
the thin films, which may suffer from the effects of sub- 
strate induced strain and/or defects present in the films. 



Moreover, substrate contribution to the magnetization 
and the weight of the material in the thin film samples 
are difficult to measure accurately. All these things may 
lead to different magnetization value than bulk sample. 
Interestingly, Matsuno et al. have explained their satura- 
tion magnetization value by assuming Co 4+ in intermedi- 
ate spin state and that would lead to the same magnetic 
moment we obtained. 



IV. CONCLUSIONS 

The detailed electronic and magnetic properties of a 
single-layered compound - Sr2Co04 - have been studied 
using ab initio calculations with different formulations of 
the density functional theory. The GGA calculation gives 
a ferromagnetic metallic state with total magnetic mo- 
ment of ~ 2 jUs/fii. The GGA+{7 calculations provide 
two ferromagnetic solutions with S = 3/2 spin config- 
uration: (i) half-metallic and (ii) metallic. The ferro- 
magnetic half-metallic state is the lowest energy state 
and representative of the ground state. Our results also 
suggest that the compound may show metal to insulator 
transition when a small amount hole is doped at the Co 
sites. In the present work we have also estimated the 
crystal-field energy, magnetic moments, Heisenberg ex- 
change parameters, transition temperature, etc., which 
provide good description of experimentally observed elec- 
tronic and magnetic properties of the system. 
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